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1.  Project  objectives. 


The  following  objectives  were  originally  proposed,  and  the  course  of  progress  during  this 
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(a)  Working  with  the  MU  radar  group  of  Kyoto  University,  describe  the  salient  features  of 
variations  of  Cn^  and  e  at  the  MU  and  WS  VHP  radars  with  height,  time  of  day,  season, 
background  winds  and  other  weather  conditions.  [11(e),  III(a.7)] 

(b)  Compare  the  intensities  of  gravity  wave  activity  at  MU  and  WS  and  their  relationship  to 
variables  of  the  large-scale  flow  such  as  jet  stream  intensity,  wind  shear,  season,  and  time-of- 
day.  [III(a.6)] 

(c)  Compare  the  intensity  of  gravity  wave  activity  from  radar  measurements  with  those  from 
other  measurement  platforms,  such  as  GPS/MET  satellite,  to  generate  an  algorithm  for 
extrapolating  results  seen  at  radar  sites  to  remote  regions  of  the  world.  [11(a)] 

(d)  Develop  conceptual,  physical,  or  statistical  models  of  the  variations  of  Cn^  or  s  as  functions 
of  atmospheric  variables  of  the  large-scale  flow  that  are  routinely  available  such  as  static 
stability,  jet  stream  features  (horizontal  proximity,  vertical  wind  shear,  maximum  speed, 
ageostrophic  wind  speed  indicators,  jet  height,  tropopause  discontinuity,  among  others),  season 
and  time  of  day.  [11(h)] 

(e)  Using  the  above  models,  test  their  effectiveness  by  using  predictions  of  the  large-scale  flow 
from  routinely  available  numerical  weather  prediction  models  to  predict  levels  of  Cn^  or  s  at  MU 
or  WS. 

(f)  Make  case  studies  of  the  intensity  of  observed  and  derived  turbulence  variables  during  these 
periods  to  draw  conclusions  regarding  relevant  physical  processes  and  relationships  to  mean 
winds,  stability,  and  background  weather  conditions.  [II(i),  III(a.l4)] 
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(a)  Nastrom,  G.D.,  A.  R.  Hansen,  T.  Tsuda,  M.  Nishida,  and  R.  Ware,  2000:  A  comparison  of 
gravity  wave  energy  observed  by  VHP  radar  and  GPS/MET  over  central  North  America.  J. 
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radar  observations  at  MU  and  White  Sands,  Ann.  Geophvs..  19,  883-888. 
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(h)  Nastrom,  G.D.,  and  F.  D.  Eaton,  2002:  A  case  study  of  wave  and  turbulence  changes  from 
MST  radar  observations  at  Vandenburg  AFB  during  cyclone  passage.  J.  Appl.  Meteor.,  in 
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(3)  Hansen,  A.  R.,  G.D.  Nastrom,  Eaton,  F.D,  1999:  The  annual  cycle  of  gravity  wave  activity  at 
5-20  km  from  VHF  radar  observations  at  White  Sands  Missile  Range.  T^ew  Mexico,  Preprints, 
29th  Conference  on  Radar  Meteorology,  AMS,  Montreal,  12-16  July,  509-512. 

(4)  Nastrom,  G.D.,  1999:  Measurements  of  eddy  dissipation  rate  and  irmer  and  outer  scales  of 
turbulence  from  MST  radar.  Invited,  NCAR  Workshop  on  meter-scale  turbulence,  Boulder, 
August  9-11. 

(5)  Hansen,  A.  R.,  G.D.  Nastrom,  Eaton,  F.D,  2000:  Seasonal  variation  of  gravity  wave  activity 
at  5  -  20  km  observed  with  the  VHF  radar  at  White  Sands  Missile  Range,  New  Mexico. 
Proceedings,  9*^  International  Workshop  on  Technical  and  Scientific  Aspects  of  MST  Radar, 
Toulouse,  March  13-17. 

(6)  Nastrom,  G.D.,  and  T.Tsuda,  2000:  Comparisons  of  turbulence  intensity  and  its  variability 
over  New  Mexico  and  Japan,  Proceedings,  9*  International  Workshop  on  Technical  and 
Scientific  Aspects  of  MST  Radar,  Toulouse,  March  13-17. 

(7)  Nastrom,  G.D.,  and  T.Tsuda,  2000:  Anisotropy  of  Doppler  spectral  parameters  in  the  MU  and 
White  Sands  VHF  radar  observations.  Proceedings,  9*'’  International  Workshop  on  Technical 
and  Scientific  Aspects  of  MST  Radar,  Toulouse,  March  13-17. 
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(8)  Nastrom,  G.D.,  J.  Furumoto,  T.E.  VanZandt,  and  T.  Tsuda,  2000:  A  Dual-Beamwidth 
Method  for  observing  atmospheric  turbulence  intensity  with  radar,  Fifth  International 
Symposium  on  Tropospheric  Profiling:  Needs  and  Technologies.  Adelaide,  December  4-8. 

(9)  Nastrom,  G.D.,  2000:  Comparisons  of  turbulence  over  Japan  and  New  Mexico  using  radar 
observations,  ABL-AFOSR  Atmospheric  Workshop,  Kirtland  AFB,  July  26-28. 

(10)  VanZandt,  T.E.,  W.  L.  Clark,  J.  Furumoto,  T.  Tsuda,  and  G.  D.  Nastrom,  2001 :  A  Dual- 
Beamwidth  Radar  Method  for  Measuring  the  Turbulent  Velocity 

Variance,  URSI  2001,  Boulder,  January. 

(11)  Nastrom,  G.D.,  J.  Furumoto,  T.  Tsuda,  T.E.  VanZandt,  C.  Love,  and  W.L.Clark,  2001 :  A 
Dual-Beamwidth  Method  for  observing  atmospheric  turbulence  intensity  with  radar,  lAMAS 
2001,  Innsbruck,  July  10-18. 

(12)  Nastrom,  G.D.,  A.  R.  Hansen,  J.  Otkin,  and  F.D.  Eaton,  2001:  MST  radar  observations  of 
gravity  waves  and  turbulence  near  thunderstorms.  LAMAS  2001,  Innsbruck,  July  10-18. 

(13)  Nastrom,  G.D.,  F.  D.  Eaton,  and  A.  R.  Hansen,  2002:  Some  sources  of  gravity  waves  as 
seen  from  MST  radar  observations.  International  Symposium  on  Equatorial  Processes  Including 
Coupling,  Kyoto,  Japan,  March  1 8-22,  Submitted. 

(14)  Nastrom,  G.  D.,  F.  D.  Eaton,  and  E.  Boll,  2002:  A  case  study  of  atmospheric  conditions  at 
2-19  km  over  Vandenburg  AFB  during  passage  of  a  cyclone.  10*"  Conference  on  Aviation, 
Range,  and  Aerospace  Meteorology,  Portland,  May  13-16,  Submitted. 

(b)  Consultative  and  advisory  functions 

At  the  request  of  Dr.  E.  M.  Dewan  of  AFRLA^SBM,  Hanscom  AFB,  MA,  informal  review 
comments  were  given  on  a  journal  article  on  gravity  waves  that  he  is  preparing.  Also,  graphs  of 
data  from  another  project  on  gravity  waves  and  turbulence  were  provided  to  Dr.  Dewan  at  his 
request..  Also  at  the  request  of  Dr.  E.  M.  Dewan  of  AFRLWSBM,  Hanscom  AFB,  MA,  a  brief 
effort  was  made  to  evaluate  a  prediction  of  his  gravity  wave  model.  As  this  task  will  require 
considerable  time  it  will  be  held  until  further  time  becomes  available. 

IV.  Other  results 

(a)  Cn  and  synoptic-scale  vorticity. 

See  Appendix  A. 

(b)  Results  from  the  Vandenburg  AFB  MST  radar. 

See  Appendix  B. 

(c)  Case  studies  of  waves  and  turbulence  at  the  MU  MST  radar. 


See  Appendix  C. 
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V.  Summary  of  accomplishments/new  findings. 

(a)  The  demonstration  that  the  gravity  wave  intensity  seen  by  satellites  (GPS/MET) 
agrees  with  that  seen  by  the  radar  at  White  Sands.  This  agreement  establishes  confidence  in  the 
GPS/MET  values;  since  the  satellite  data  are  global  this  result  provides  a  means  of  studying  the 
global  distribution  of  gravity  wave  intensity. 

(b)  The  demonstration  that  gravity  wave  source  mechanisms  are  different  at  WS  and  MU 
although  the  wave  energies  are  about  the  same  magnitude. 

(c)  The  demonstration  that  Cn^  in  the  stratosphere,  above  18  km,  is  not  changed  by  the 
presence  of  thunderstorms  in  the  troposphere. 

(d)  The  demonstration  that  the  diurnal  changes  in  Cn^  and  spectral  width  above  the  mid¬ 
troposphere  are  very  small  during  all  seasons  at  both  WS  and  MU. 

(e)  The  demonstration  that  Cn^  in  the  stratosphere  at  both  WS  and  MU  is  very  highly 
correlated  with  relative  vorticity.  The  correlation  is  improved  when  distance  from  the 
tropopause  is  used  as  the  vertical  coordinate  and  when  the  is  normalized  by  local 
atmospheric  density.  The  model  developed  is  thus  a  combined  conceptual/  physical/statistical 
model. 

(f)  The  demonstration  in  the  case  study  at  Vandenburg  AFB  that  turbulence  (Cn^  and 
specific  turbulent  kinetic  energy)  and  waves  increase  during  the  passage  of  a  cyclone  are 
generated  only  by  upper  tropospheric  effects  because  upwelling  effects  from  the  planetary 
boundary  layer  are  cut  off  by  a  strong  surface  inversion. 

(g)  The  demonstration  in  the  case  studies  at  MU  that  turbulence  intensities  in  the  lower 
stratosphere  are  often  more  closely  related  to  the  intensity  of  gravity  waves  than  to  in  situ  wind 
speed  and  large-scale  wind  vertical  shear. 

VI.  Personnel. 

The  following  persons  were  supported  during  this  effort. 

Dr.  Gregory  D.  Nastrom,  PI 

Dr.  Anthony  R.  Hansen,  (Augsburg  College),  Contractor 


Appendix  A.  Cn  and  the  Synoptic-Scale  Vorticity 

To  determine  if  a  linkage  exists  between  the  intensity  of  mesoscale  turbulence  and 
synoptic-scale  features  in  the  mid  latitude  atmosphere,  a  study  was  undertaken  to 
compare  radar  turbulence  measurements  with  analyzed  synoptic  scale  wind  fields.  The 
radar  data  used  in  this  study  consisted  of  the  uncalibrated  back-scattered  power  measured 
at  the  MU  radar  in  Japan  from  1991-1996  and  Cn^  measurements  taken  at  the  WSMR 
radar  during  the'same  period.  The  analysis  period  was  confined  to  the  cold  season 
months  when  mean  wind  speeds  in  the  upper  troposphere  were  relatively  high  at  both 
stations.  Months  used  for  study  at  MU  were  October  through  April  of  each  year  and  at 
WSMR  the  months  used  were  November  through  April  of  each  year.  Analysis  of 
synoptic-scale  meteorological  fields  was  accomplished  through  the  use  of  NCAR/NCEP 
reanalysis  data  for  the  same  periods.  This  report  will  concentrate  on  comparisons 
involving  analyzed  winds  at  300  mb  and  500  mb.  The  analyzed  winds  were  used  to 
compute  wind  shears  and  relative  vorticity  of  the  synoptic  scale  flow.  The  objective  is  to 
establish  what  relationship  exists  between  variations  in  the  synoptic-scale  winds  and  the 
locally  measured  turbulence  in  the  vicinity  of  the  tropopause. 

Consider  first  the  results  obtained  using  the  MU  data.  Winds  for  the  grid  points 
surroimding  the  MU  radar  site  were  extracted  from  the  NCAR/NCEP  reanalysis  and  the 
relative  vorticity,  for  the  location  directly  above  the  radar  was  computed.  The 
reanalysis  data  were  available  at  6  hour  intervals.  The  resulting  relative  vorticity  values 
were  compared  to  the  hourly  means  of  the  uncalibrated,  back-scattered  power 
measurements  from  MU.  (The  uncalibrated,  back-scattered  power  is  used  as  a  proxy  for 
Cn^.)  Vorticity  values  were  considered  constant  within  3  hours  of  a  reanalysis  time  for 
the  purposes  of  computing  correlations  with  the  hourly  back-scattered  power 
measurements.  In  addition,  to  remove  the  possibility  that  any  resulting  correlation  was 
the  result  merely  of  vertical  displacements  of  the  tropopause  resulting  from  the  passage 
of  synoptic-scale  waves  or  jet  streaks,  the  measured  back-scattered  power  profiles  were 
adjusted  vertically  to  a  tropopause  relative  co-ordinate  system.  To  establish  the 
tropopause  height  above  MU,  routine  radiosonde  measurements  from  the  nearby  stations 
of  Wajima  (37.22N  136.54E)  and  Shiono-misaki  (33.27N,  135.46E)  were  used. 
Tropopause  heights  based  upon  the  WMO  definition  were  computed  from  the  soundings 
at  12  hourly  intervals  and  tropopause  heights  for  MU  were  interpolated  in  space  and  time 
from  these  results.  Thus,  for  each  hourly  mean  radar  measurement,  a  tropopause  height 
was  determined  and  the  radar  power  profile  was  adjusted,  based  up  the  mean  profile  from 
the  dataset,  to  a  tropopause  relative  vertical  co-ordinate.  These  values  were  then 
correlated  with  the  MU  vorticity  values  determined  from  the  reanalysis. 

In  addition,  a  factor  to  normalize  the  radar  data  to  remove  possible  dependence  of 
the  backscatter  power  values  on  the  air  density  was  employed.  Measurements  of  the  air 
pressure  and  temperature  at  the  level  of  the  WMO-definition  tropopause  for  MU  were 
computed  from  the  radiosonde  data.  From  these  a  factor  to  normalize  these  results  from 
the  air  density  at  the  tropopause  for  each  hour  to  the  air  density  at  14  km  was  employed. 
This  factor 

N  —  (ptrop/pi4km) 


is  designed  to  remove  any  possible  dependence  in  the  results  on  mere  density  differences 
between  the  tropopause  level  on  one  day  versus  another.  A  selection  of  results  is 
presented  below. 

Fig.  A-1  illustrates  the  comparison  of  all  hoxirly  power  and  vorticity  pairs  for  MU, 
using  the  relative  vorticity  at  300  mb  at  an  altitude  1 .2  kin  below  the  local  tropopause. 

The  relationship  at  this  altitude  was  relatively  strong.  At  higher  altitudes,  particularly 
above  the  tropopause,  similar  results  also  appeared  in  that  the  regression  curves  between 
the  relative  vorticity  and  the  departure  of  the  back-scattered  power  from  its  mean  are 
nearly  identical,  but  the  correlations  between  the  regression  line  fit  to  the  data  and  the 
data  is  lower  at  altitudes  above  the  tropopause.  Inspection  of  this  figure  clearly  indicates 
a  relationship  between  increasing  cyclonic  vorticity  and  increasing  back-scattered  power 
(and  therefore  increased  levels  of  turbulence).  The  correlation  of  r=0.41  between  this 
data  sample  and  the  plotted  linear  regression  fit  to  the  data  is  quite  impressive 
considering  the  noise  level  in  the  underlying  data.  To  reduce  Ae  noise  level  in  this  plot, 
the  data  pairs  were  binned  within  ranges  of  relative  vorticity  and  means  computed.  These 
results  are  presented  in  Fig.  A-2.  Error  bars  in  this  figure  correspond  to  estimates  of  the 
standard  error  in  the  mean  of  the  data  plotted.  The  resulting  correlation  of  a  linear 
regression  fit  to  the  averaged  data  is  quite  high. 

Corresponding  results  at  500  mb  for  MU  are  similar  to  those  at  300  mb  but  fit  of 
the  regression  lines  to  the  data  are  not  as  strong.  Fig.  A-3  illustrates  the  same  results  as 
Fig.  A-2  except  that  the  vorticity  at  500  mb  was  used. 

The  same  procedure  was  followed  for  the  WSMR  radar  data.  Reanalysis  data 
corresponding  to  the  location  of  WSMR  were  extracted  and  vorticity  values  computed. 
The  measured  Cn^  values  from  the  WSMR  radar  were  also  adjusted  vertically  to  produce 
tropopause  relative  co-ordinates.  The  WSMR  tropopause  was  determined  from 
radiosonde  measurements  from  nearby  El  Paso,  TX  (ELP).  In  addition,  the  same  density 
normalizing  factor,  N,  was  used  to  remove  any  dependence  of  the  results  on  mere  altitude 
differences  between  the  tropopause-relative  vertical  co-ordinate  at  one  hour  versus  that  at 
another  hour. 

The  results  for  WSMR,  also  at  1.2  km  below  the  tropopause  between  the 
normalized  Cn^  and  the  relative  vorticity  at  500  mb  is  given  in  Fig.  A-4.  A  very  similar 
relationship  to  that  found  at  MU  is  also  apparent  in  this  independent  data  sample.  That  is, 
normalized  Cn^  values  increase  with  increasingly  cyclonic  relative  vorticity  in  the 
synoptic-scale  flow.  For  anticyclonic  conditions,  the  Cn^  values  are  relatively  lower.  For 
cyclonic  conditions,  Cn^  values  are  higher.  Results  based  purely  on  the  lateral  shear  of 
the  synoptic  scale  flow  lead  to  essentially  the  same  results  (not  shown).  The  results  at 
WSMR  show  a  better  fit  to  the  regression  curve  at  500  mb  than  do  the  corresponding 
results  at  MU. 
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Appendix  B.  Results  from  the  Vandenburg  AFB  MST  radar 

A  case  study  of  atmospheric  conditions  at  2-19  km  over  Vandenburg 
AFB  during  passage  of  a  cyclone 


I.  Introduction. 

The  purpose  of  this  report  is  to  describe  a  case  study  of  winds  and  turbulence  over 
Vandenburg  AFB,  CA,  during  7/20/20Z  to  to/23Z  July  2001.  During  this  period  *e 
50  MHz  radar  at  Vandenburg  AFB  was  operated  continuously.  Vertical  profiles  of  the 
wind  speed,  backscattered  power  (reported  as  signal-to-noise  ratio,  SNR,  and  calibrated 
as  the  refractivity  turbulence  structure  constant,  Q,^),  and  the  Doppler  spectral  width 
(HPHW  in  ms’’)  were  obtained  at  3-niinute  intervals  on  the  three  beams  (vertical  and  1 5 
degrees  from  the  zenith  at  azimuths  225  and  315  degrees).  The  VBG  radar  is  similar  m 
size  (half-power  beamwidth  3.2  degrees)  and  power  (208  kW  peak)  to  that  operated  at 
White  Sands  Missile  Range  during  the  1990s.  The  radar  data  are  augmented  by 
rawinsonde  ascents  and  surface  weather  observations  at  Vandenburg  AFB  and  by  routine 
weather  map  analyses. 

II.  Data  Quality 

A.  Initial  data  inspection.  At  first  glance  the  vertical  profiles  of  the  observations 
look  reasonable.  For  example,  Fig.  1  shows  the  individual  profiles  of  log  Q, ,  SNR  (in 
dB),  velocity  (ms'*),  and  log  spectral  width  for  each  beam  (beams  vertical,  315  ,  and  225 
°  in  rows  upper,  center,  and  lower,  respectively)  for  the  hour  7/21/OOZ.  The  number  of 
profiles  available  (16  or  20)  is  given  above  the  Cn^  plot  for  each  beam.  A  vigorous 
quality  control  algorithm  was  used;  briefly,  Cn^,  velocity,  and  width  were  not  used  if 
SNR<-12  dB  or  if  velocity  or  failed  temporal  consistency  checks  or  if  the  vertical 
shear  of  horizontal  wind  exceeded  realistic  limits.  Note  that  the  apparent  layered 
structure  of  is  persistent  throughout  the  hour.  Also,  many  of  Ae  small-scale  features 
of  the  horizontal  velocity  profiles  are  persistent  (horizontal  velocities  were  estimated  by 
multiplying  the  radial  velocities  by  esc  15°;  i.e.,  the  effects  of  vertical  velocities  were 
assumed  to  be  negligible).  Although  the  relative  variability  of  the  width  profiles  at  a 
given  height  is  larger  than  that  of  Cn^  and  wind,  the  general  shape  appears  constant 
throughout  the  hour.  Inspection  of  corresponding  plots  from  the  odier  hours  (not  shown) 
gives  similar  results.  Thus,  the  data  appear  reasonable;  for  convenience,  hourly  medians 
of  the  individual  observations  at  each  level  will  be  used  below.  Any  minor  adjustments, 
such  as  re-calibration  of  the  Cn^  values,  will  not  change  the  present  conclusions. 

B.  Sensitivity  limits.  Note  in  Fig.  1  that  the  left-hand  side  of  the  oblique-beam 
Cn^  profiles  appear  to  be  “clipped”,  that  is,  there  seems  to  be  a  lower  limit  of  sensitivity 
that  slopes  upward  toward  the  right.  This  feature  is  confirmed  in  Fig.  2.  Fig.  2  shows  Ae 
hourly-median  profiles  for  all  172  hours  in  the  lower  panels.  The  upper  panels  show  the 
profile  of  the  median  of  the  172  hourly  values  (solid)  and  the  profiles  of  the  maximum 
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and  minimum  value  at  each  level.  Also  plotted  in  Fig.  2  are  the  sensitivity  IMts  (called 
the  “minimum  detectable"  limit  by  Green  et  al.,  1979).  The  nummum 
computed  by  multiplying  the  reported  by  the  ratio  SN^in/SNRr,^rted,  where  SN^,n 
was  taken  to  be  -12  dB.  The  “clipped”  left  edge  of  the  Cn^  profiles  obviously  anses  from 
the  effects  of  sensitivity  limits.  Also,  note  that  the  minimum  detectable  signal  curve  has 
a  knee  of  uncertain  origin  at  about  3.5  km;  thus,  we  should  not  use  Cn  values  below  3.5 

km.  ■ 


C.  Anisotropy.  Fig.  3  compares  the  median  profiles  of  G„^  and  Width  from  the 

oblique  beams.  Right  hand  panels  show  the  ratios  of  Beam  225  to  Bearn  315.  The  vdues 
of  Cn'  from  Beam  225  are  larger  (from  0  to  2  dB)  in  the  troposphere,  and  above  about  14 
km  the  values  are  indistinguishable.  In  view  of  the  very  small  differences  at  high 
altitude,  where  tilting  of  specular  layers  might  be  quite  small,  this  result  suggests  that  t  e 
two  beams  are  well-calibrated  with  respect  to  each  other.  Regarding  width,  the  values  for 
Beam  315  are  very  slightly  larger  above  4  km.  However,  note  that  the  plotted  values  are 
observed  widths  and  that  corrections  for  beam-,  shear-,  and  wave-broademng  ettects 
must  be  applied  to  these  (e.g.,  Nastrom  and  Eaton,  1997).  In  summary,  the  two  oblique 
beams  appear  to  give  similar  results. 


III.  Radar  Observations 

A.  Radar  wind  observations.  Fig.  4  (upper  panels)  shows  the  hourly  meians  of 
zonal  and  meridional  wind  speeds  (ms'*).  For  reference,  hour  100  is  7/25/OOZ.  color 
contours  are  at  10  ms’*  intervals;  for  zonal  wind  the  zero  line  is  that  betw^  light  blue 
and  dark  blue  and  for  meridional  wind  it  is  that  between  blue  and  green.  The  salient 
features  of  zonal  wind  are  the  westward  winds  below  about  6  km  from  hours  90-1^50,  the 
waxing  and  waning  of  the  jet  near  13  km  with  strongest  pulse  near  hour  90,  and  the 
quasi-periodic  oscillations  above  16  km  with  downward  phase  progression  with  tune. 

The  main  features  of  meridional  wind  are  the  strong  burst  of  scatbward  winds  near  9  Km 
at  hour  100  and  the  quasi-periodic  oscillation  above  16  km.  Conespondmg  wmd  speed 
and  wind  direction  charts  are  given  in  the  lower  panels  where  each  of  the  features 

mentioned  above  can  easily  be  found.  -n  u 

The  quasi-periodic  wind  behavior  above  16  km  is  interestog;  however,  it  will  be^ 

seen  in  Figs.  5  and  6  below  that  it  does  not  seem  to  have  a  large  influence  on  SNR  or  C„  . 
Since  this  report  concerns  Cn',  the  periodic  oscillation  will  not  be  discussed  further. 

B.  Radar  SNR.  Fig.  5  shows  the  time-height  plots  of  SNR  in  dB.  Those  for 
Beam  315  (225)  are  in  the  upper  (lower)  panel.  The  left-hand  panel  shows  the  c^ttire 
range  2-21  km.  Very  large  SNR  values  are  seen  in  the  troposphere,  below  about  6  km. 
At  higher  altitudes  the  main  feature  is  an  enhancement  of  SNR  at  about  8-16  km  from 
about  hour  70-1 12  (about  7/23/1 8Z  to  7/25/12Z).  This  enhancement  is  seen  m  more 
detail  in  the  right-hand  panels  which  show  only  the  region  above  10  km. 

C.  Radar  Cn'.  Fig.  6  is  similar  to  Fig.  5,  except  it  shows  log  Cn'.  The  logarithm 
is  used  hereafter,  but  will  be  referred  to  as  just  Cn's  because  past  studies  have  shown  the 
logarithm  of  Cn'  is  normally  distributed.  The  main  features  in  Fig.  6  are  the  very  rapid 
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decrease  of  Cn^  above  about  6  km  and  the  apparent  enhancement  of  Cn  up  to  about  16 
km  during  approximately  hours  70-1 12.  Due  to  the  large  range  of  Cn  values  m  the  ett- 
hand  panels  of  Fig.  6  it  is  difficult  to  tell  how  low  this  feature  extends^  To  help  see  the 
extent  of  the  enhancement,  Fig.  7  shows  plots  of  the  perturbations  of  Cn  ,  defined  as  the 

hourly  values  minus  the  median  of  all  hourly  values  at  each  level.  Profiles  of  die 

medics  are  given  in  the  upper  panel  for  convenience  (Beam  315  is  the 

that  the  enhancement  in  Fig.  7  (especially  for  Beam  225)  exten^  from  about  hour  70-112 

(about  7/23/1 8Z  to  7/25/1 2Z)  and  that  it  extends  from  below  4  km  (the  mimmum  level  ot 

reliable  Cn^  values)  to  about  15  km. 

D.  Turbulent  Kinkic  Energy  (TKE).  TheupperpanelsofFjg  8  show  the  radar 

spectral  width  plotted  as  the  turbulent  kinetic  energy  density,  o,  (m  s  ).  These  values 
have  been  corrected  for  beam-,  shear-,  and  wave-broademng  effects  followmg  Nasfrom 
and  Eaton  (1997).  The  main  feature  is  the  general  enhancement  of  o,  from  about  hour 
70-1 12.  This  time  period  correlates  very  well  with  the  enhancement  times  noted  for  SN 
and  Cn^  although  the  relative  intensity  of  the  enhancement  of  Ot  is  smaller  than  that  for 

SNR  or  Cnl 

E.  Standard  deviation  of  vertical  velocity.  The  lower-right  panel  of  Fig.  8  shows 
the  hourly  standard  deviations  of  vertical  velocity,  Ow  (ms' ).  The  values  ^e  genera  y 
very  small,  less  than  about  0.1  ms'*,  except  between  about  hours  80  and  1 10  at  about  9-14 
km  when  values  exceed  0.25  ms  . 

rV.  Case  study  results 

A.  Synoptic  weather  patterns.  Fig.  9  shows  the  500  mb  maps  at  12Z  daily  for 
7/20-7/27  2001.  The  flow  over  Vanderiburg  AFB  was  from  the  southwest  through  7/22. 
On  7/23  a  sharp  trough  began  to  influence  conditions  over  southern  California  as  win^  at 
VBG  became  westerly.  The  trough  deepened  and  formed  into  a  closed  low  centered  just 
north  of  VBG  on  7/24.  By  7/25  the  low  had  filled  a  little  and  its  center  had  moved 
several  hundred  km  east  of  VBG.  By  7/26  the  remaining  weak  trough  was  over  Arizona 
with  the  flow  over  VBG  much  weaker.  Flow  was  back  to  southwest  on  7/27. 

Rawinsonde  observations  from  VBG  for  23/00Z-27/12Z  are  given  m  Fig- 10 
launch  was  made  from  South  Vandenburg  on  25/OOZ,  a  site  about  10  km  soufo  of  yBG). 
The  temperature  soundings  show  a  strong  inversion  over  VBG  throughout  this  penw 
with  the  top  of  the  inversion  several  thousand  feet  above  the  surface.  Surface  wither 
observations  from  VBG  (not  shown)  indicate  continuous  low  cloud  coverage  wifli 
occasional  fog  throughout  this  period.  All  local  hills  were  well  below  the  top  of  the 
inversion,  and  winds  in  the  lower  levels  were  consistently  relatively  hght,  10  or 
less,  throughout  this  period.  Thus,  it  appears  that  the  upper  levels  were  isolated  from  the 

boundary  layer.  u- 

Note  that  the  winds  from  the  radar  (Fig.  3)  agree  very  well  with  the  geostrophic 

flow  pattern  (Fig.  9)  and  the  winds  from  the  rawinsondes  (Fig.  10).  This  agreement  gives 

further  confidence  in  the  radar  observations. 

At  first  glance,  the  strong  wind  changes  in  the  upper  troposphere  seen  in  Fig.  4,  as 
well  as  the  enhancements  in  SNR,  Cn^,  and  Ow  in  Figs.  5-8,  will  be  attribut  to  e 
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passage  of  the  upper  level  cyclonic  circulation  feature.  Due  to  the  strong  low  leve 
inversion  and  light  low  level  winds,  it  seems  unlikely  that  boundary  layer  effects  (such  as 
flow  over  rough  terrain  or  low  level  convection)  were  the  source  of  the  enhancements  in 
SNR,  Cn^  and  Ow  seen  in  Figs.  5-8. 

B.  Correlation  patterns.  Fig.  11  shows  time  series  of  the  median  Cn^Ot^,  and  Ow 

over  all  levels  in  the  layer  10-15  km  in  the  left-hand  panels.  The  right-hand  panels  show 
plots  of  the  same  data,  scaled  to  have  a  range  of  unity.  The  upper  (center)  p^e 
compares  Cn^  with  (Ow)  and  the  lower  panel  compares  Ot  and  Ow  There  is  high 
correlation  of  Cn^  with  (0.87).  The  correlation  of  C„^  with  Ow  is  also  fairly  high 
(0.60),  and  note  that  o,^  and  Ow  are  fairly  well-correlated  (0.59). 

The  high  correlation  coefficient  for  C„'  with  o.'  (0.87)  is  important  because  it 
implies  that  a  large  part  of  the  variability  in  Cn^  is  due  to  turbulence  changes.  In  genera  , 
local  values  of  depend  on  both  the  intensity  of  turbulence  through  the  outer  scale  of 
turbulence,  U,  and  the  mean  gradient  of  refractive  index  (M) 

Cl=a^a'LfM^  W 

where  a^  is  a  constant  and  a’  is  the  ratio  of  eddy  diffusivities,  usually  taken  as  constant 
(see  Gage,  1990,  for  a  review).  Lo  is  is  related  to  the  eddy  dissipation  rate,  £,  and  the  win 
shear,  du/dz,  by  Lo=(e/Y)’^ldu/dzr^'^  where  y  is  usually  taken  as  constant.  M  is  related  to 
the  stability  (dlnO/dz)  and  specific  humidity  (q)  as 


M=-77xl0'‘^ 

T 


d\n0  15500g  dine  7800  dq\ 
dz  T  dz  T  dz) 


(2) 


In  the  presence  of  high  humidity  and  humidity  gradients  M  and  hence  Cn  depends 
strongly  upon  humidity  (Tsudaetal.,  1988).  However,  the  rawinsonde  reports  durmg 
this  period  (Fig.  10)  show  that  humidity  levels  were  low  and  uniform  with  time  above 
about  700  mb  (about  3  km).  The  rawinsonde  reports  further  show  that  the  stability  was 
uniform  with  time.  Nastrom  and  Eaton  (2001)  show  that  in  the  absence  of  humidity  the 
volume  averaged  is  proportional  only  to  turbulence  intensity  and  stability,  and  since 
the  stability  changes  here  were  small  the  changes  in  Cn^  with  time  were  due  to  changes  m 
turbulence  intensity.  The  high  correlation  of  with  Ot  seen  in  Fig.  1 1  supports  this 
suggestion. 

2 

C.  Source  of  turbulence.  Fig.  12  compares  time  series  of  the  hourly  median  C„ 
at  10-15  km  with  those  of  wind  speed  and  vertical  wind  shear  over  300  m  intervals.  (For 
reference,  time  series  of  the  median  Cn^  and  median  wind  speed  and  direction  over  1  km 
intervals  at  10-15  km  are  given  in  Fig.  13.  Note  in  Fig.  13  tiiat  the  wind  speed  and 
direction  are  very  similar  at  all  levels  10-13  km  and  are  different  at  14-15  km.) 

Returning  to  Fig.  12,  note  that  the  correlation  between  Cn  and  wind  shear  (0.80)  is  much 
greater  than  that  between  Cn^  and  wind  speed  (0.49),  suggesting  that  the  turbulence  is 
generated  by  a  local  dynamic  instability  and  that  this  instability  is  not  necessanly  related 
to  wind  speed.  That  is,  since  Cn^  is  not  closely  related  to  wind  speed  (wind  speed  is  an 
indicator  of  the  large-scale  flow),  the  turbulence  may  be  due  to  instabilities  in  small-scale 
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flow  features  such  as  gravity  waves  (this  suggestion  is  supported  by  the  correlation  of  Cn 
and  Ow  in  Fig.  1 1  since  Ow  is  generally  regarded  as  an  indicator  of  gravity  wave  intensity). 

Possible  sources  of  gravity  waves  include  geostrophic  adjustment  processes  and 
shear  instabilities  along  the  edges  of  the  jet  stream.  Since  the  flow  is  cyclonic  it  seems 
unlikely  that  inertial  instability  processes  are  important,  and  since  the  planetary  boundary 
layer  is  capped  by  a  strong  inversion  with  weak  low-level  flow  it  seems  unlikely  that 
upwelling  of  waves  from  the  low  levels  is  important. 

In  another  Appendix  of  this  report  it  is  shown  that  Cn^  at  White  Sands  Missile 
Range  and  retained  power  at  MU  in  Japan  are  well-correlated  with  vorticity.  Fig.  14 
shows  the  time  series  of  geostrophic  vorticity  at  200  mb  over  VBG  computed  using  the 
height  field  shown  on  routine  weather  maps.  The  main  feature  is  a  spike  in  vorticity 
centered  at  25/OOZ;  this  time  corresponds  to  hour  100  in  Fi^.  4-7  and  thus  is  about  at  the 
center  of  the  period  of  enhanced  Cn  .  This  case  study  provides  further  evidence  for  the 
usefulness  of  vorticity  as  a  predictor  of  Cn^. 
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SLAT  34.75 
SLON  -120. 
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SHOW  12.10 
LIFT  15.64 
LFTV  15.60 
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SHOW  4.94 
LIFT  11.23 
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SLAT  3475 
SLON  -120. 
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SHOW  2.36 
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PWAT  25.23 


10 


72393  VBG  Vandenberg  Afb 


SLAT  34.75 
SLON  -120. 
SELV  121.0 
SHOW  6.15 
LIFT  12.70 
LFTV  12.58 
SWET  61.48 
KJNX  11.10 
CTOT  7.70 
VTOT  30.70 
TOTL  38.40 
CAPE  0.00 
CAPV  0.00 
CINS  0.00 
CINV  0.00 
EQLV  -9999 
EQTV  -9999 
LFCT  -9999 
LFCV  -9999 
BRCH  0.00 
BRCV  0.00 
LCLT  285.0 
LCLP  950.4 
MLTH  289.2 
MLMR  9.32 
TUCK  5786. 
PWAT  22.07 


OOZ  26  Jul  2001 


72393  VBG  Vandenberg  Afb 


SLAT  34.75 
SLON  -120. 


SELV  121.0 
SHOW  5.24 
LIFT  15.07 
LFTV  15.01 
SWET  67.54 
KJNX  15.70 
CTOT  9.50 
VTOT  30.50 
TOTL  40.00 
CAPE  0.00 
CAPV  0.00 
CINS  0.00 
CINV  0.00 
EQLV  951.6 
EQTV  951.7 
LFCT  953.7 
LFCV  953.7 
BRCH  0.00 
BRCV  0.00 
LCLT  283.7 
LCLP  953.7 
MLTH  267.6 
MLMR  8.51 
THCK  5777. 
PWAT  24.35 


-40  -30  -20  -10  0  10  20  30  40 


12Z  26  Jul  2001 


10 


72393  VBG  Vandenberg  Afb 
100 


200 


300 


400 

500 

600 

700 

800 

900 


\ili_ 


'A 


f 


SLAT 

34.75 

SLON 

-120. 

SELV 

121.0 

SHOW  2.93 

LIFT 

12.47 

LFTV 

12.41 

SWET 

106.8 

KINX 

22.10 

CTOT 

12.30 

VTOT 

31.30 

TOTL 

43.60 

CAPE 

0.00 

CAPV 

0.00 

CINS 

0.00 

CINV 

0.00 

EQLV 

-9999 

EQTV 

-9999 

LFCT 

-9999 

LFCV 

-9999 

BRCH 

0.00 

BRCV 

0.00 

LCLT 

284.5 

LCLP 

925.3 

MLTH 

290.9 

MLMR  9.31 

THCK 

5806. 

PWAT 

30.00 

OOZ  27  Jul  2001 


-40  -30  -20  -10  0  10  20  30  40 


SLAT  34.75 
SLON  -120. 
SELV  121.0 
SHOW  4.22 
LIFT  18.87 
LFTV  18.87 
SWET  86.74 
KJNX  16.10 
CTOT  10.30 
VTOT  31.30 
TOIL  41.60 
CAPE  0.00 
CAPV  0.00 
CINS  0.00 
CINV  0.00 
EQLV  -9939 
EQTV  -9999 
LFCT  -9999 
LFCV  -9999 
BRCH  0.00 
BRCV  0.00 
LCLT  279.7 
LCLP  908.9 
MLTH  287.5 
MLMR  7.36 
THCK  5796. 
PWAT  22.26 


12Z  27  Jul  2001 


10 


O' I  8  0  9  0  t-'O  30  O'O  O'L  B'O  9'0  fO  3’0  O'O  O' I  8'0  9'0  VO  Z'O  O' 
S3mVA  3MiV13y  SSmVA  3AliV13a  S3mVA  3Ali\n3y 


0'6l-  B'6l-  9'03-  6'0  L'O  B'O  C'O  t'O  H  O  03  0  91'0  310  80  0  i'O  O 

(UJ>I  gi-Ol)  P  6o|  aJVA 


48  60  72  84  96  120  144-  168  °  0  12  24  36  48  60  72  84  96  120  144  168 

VAFB  Hours  after  7/20/01  20Z  VAFB  Hours  after  7/20/01  20Z 


(UJ^I  51-01)  jO  6o|  ajVA  (|-S^0  (^>1  Ql-oi)  PSJ'ds  pu|M 


48  60  72  84  96  120  144  168  o  0  1  2  24  35  48  60  72  84  96  120  144  168 

VAFB  Hours  after  7/20/01  20Z  ’  VAFB  Hours  after  7/20/01  20Z 


61-  9'6l-  9'03- 

533^1^3  6o| 


99  93  03  31  9 

paads  puiM 


Appendix  C.  Case  studies  of  waves  and  turbulence  at  the  MU  MST  radar. 

Observations  taken  with  the  Middle  and  Upper  (MU)  Atmosphere  MST  radar  located 
near  Shigaraki,  Japan,  were  available  for  about  100  hours  per  month  during  April  1986  through 
December  1996,  as  well  as  for  an  extended  period  (422  hours  long)  during  April  1995.  The  MU 
radar  has  five  beams,  one  vertical  and  one  10  degrees  from  the  vertical  toward  each  of  the  four 
cardinal  directions.  Following  quality  control  checking,  the  individual  6-minute  profiles  were 
used  to  make  hourly  medians  of  the  backscattered  power  on  the  meridonal  beams,  horizontal 
wind  components,  and  specific  turbulent  kinetic  energy  on  the  meridional  beams,  plus  the  hourly 
standard  deviation  of  the  vertical  velocity.  Relative  backscattered  power  is  used  since  the  MU 
radar  is  not  calibrated  to  report  Cn^  and  transmitted  power  records  were  not  available  to  permit  a 
post-observation  calibration,  such  as  comparisons  with  cosmic  noise.  The  radar  data  were 
augmented  with  rawinsonde  ascents  from  Wajima  and  Shionomisaki,  located  north  and  south  of 
MU,  respectively. 

Four  cases  were  chosen  for  presentation  here,  to  illustrate  the  relationships  seen  among 
the  radar  variables. 

April  1995.  This  segment  of  observations  lasted  422  hours  and  is  the  longest  one  available. 
Time-height  plots  of  the  power  and  horizontal  wind  speed  are  given  in  Fig.  Cl .  Tropopause 
heights  at  Shionomisaki  (Wajima)  are  shown  by  black  dots  (diamonds).  The  abscissa  is  hours 
after  the  first  hour  of  the  segment.  Note  that  during  relatively  quiet  periods,  such  as  around  hour 
210  and  hour  325,  when  the  tropopause  height  shows  little  change  with  latitude,  there  is  a  sharp  , 
increase  of  power  just  above  the  tropopause  whereas  during  other  periods  the  vertical  gradient  of 
power  is  more  diffiise.  The  periods  of  large  tropopause  height  gradient  usually  correspond  with 
the  passages  of  jet  streaks,  i.e.,  pulses  of  very  strong  winds.  A  relationship  between  wind  speed 
and  power  (or  Cn^)  was  reported  by  others  (e.g.,  Hufiiagel,1974;  Gage  et  al.,  1980).  The 
individual  hourly  profiles  of  power  are  given  in  Fig.  C2,  along  with  the  median  of  all  hourly 
values  at  each  level  (crosses)  and  the  median  of  all  hourly  horizontal  wind  speeds  (dots,  shifted 
by  80  m/s  relative  to  the  values  on  the  abscissa).  Note  that  the  strongest  winds  are  at  about  12 
km,  this  level  corresponds  roughly  to  the  average  tropopause  height  (compare  Fig.  Cl).  Fig.  C4 
shows  a  time  series  of  the  mechan  of  the  hourly  values  of  power  over  all  levels  from  12  to  18  km 
(roughly,  the  lowest  6  km  of  the  stratosphere)  in  the  upper  left  panel.  Corresponding  time  series 
of  the  wind  speed  and  wind  shear  (defined  as  [(du/dz)^+(dv/dz)^]'^  where  dz=300m)  are  given  in 
the  left-center  and  lower  panels,  respectively.  Pairs  of  these  time  series  are  compared  in  the 
right-hand  panels,  where  each  is  plotted  with  a  range  of  unity.  The  correlation  coefficient  for 
each  pair  is  given  in  each  panel;  they  range  from  0.20  for  the  power-shear  pair  to  0.00  for  the 
speed-shear  pair  with  an  intermediate  value  of  0.17  for  power-speed.  Fig.  C4  is  similar  to  Fig. 

C3  except  the  variables  are  power,  specific  kinetic  energy  (c^,  determined  using  the  spectral 
width  corrected  for  beam,  shear,  and  wave-broadening  effects)  and  the  standard  deviation  of 
vertical  velocity  (Ow)  which  is  taken  as  an  indicator  of  gravity  wave  intensity.  The  largest 
correlation  in  Fig.  C4  is  that  for  the  power-Ow  pair  (0.41),  the  smallest  is  for  Ot^-Ow  (0.18),  and 
that  for  power-  (0.3 1)  is  intermediate.  The  correlation  between  power  and  Ot  indicates  that 
strong  refractivity  turbulence  is  accompanied  by  strong  mechanical  turbulence,  as  expected.  For 
this  case  study  the  closest  correlation  is  that  between  power  and  Ow,  implying  that  the  intensity  of 
high-frequency  gravity  waves  is  critical  for  the  production  of  refractivity  turbulence.  The  mean 
wind  speed  is  seen  in  this  case  to  be  a  relatively  poor  predictor  of  power,  contrary  to  HufnageTs 
model;  this  poor  positive  correlation  between  power  and  wind  speed  will  be  shown  more 
dramatically  in  two  of  the  following  three  case  studies. 

January  1996.  Figs.  C5-C8  are  similar  to  Figs.  C1-C4,  except  they  are  for  a  period  of  103  hours 
beginning  on  8  January  1996.  Also,  Fig.  C5  shows  time-height  plots  of  the  meridional  wind 
speed  (an  indicator  of  the  passage  of  troughs  and  ridges),  the  vertical  shear  per  300m,  Ow,  and  Ot . 
Fig.  C5  shows  that  the  strong  winds  early  in  the  period  were  accompanied  by  strong  Ow ,  and  the 
weak  winds  late  in  the  period  were  accompanied  by  weak  a*.  In  Fig.  C7  the  correlation  between 


power  and  wind  speed  is  seen  to  be  0.47;  while  this  value  is  large  it  is  not  as  great  as  that  for 
power-shear  (0.63),  power-  (0.59),  or  power-  Ow  (0.73). 

February  1987.  Figs.  C9-C12  are  similar  to  Figs.  C5-C8,  except  they  are  for  a  79-hour  period 
beginning  on  3  February  1987.  The  strongest  winds  are  late  in  the  period,  while  strongest  Ow  are 
early  in  the  period.  Figs.  Cl  1  and  C12  show  that  the  power  correlates  positively  with  Ow  (0.47) 
in  this  case,  where?.o.  power  porrelates  negatively  with  wind  speed  (-0.26).  The  correlation  of 
power  with  vertical  wind  shear  is  positive  but  relatively  weak  (0.17). 

May  1996.  Figs.  C13-C16  are  similar  to  Figs.  C5-C8,  except  they  are  for  a  103 -hour  period 
beginning  on  20  May  1996.  In  Fig.  C13,  wind  speed  near  12  km  shows  two  maxima  (similar  to 
Fig.  C9),  but  in  this  case  the  meridional  wind  is  first  strongly  positive,  then  negative,  then  again 
positive,  implying  the  passage  of  a  trough-ridge  pattern  (in  contrast  to  Fig.  C9  where  the 
meridional  wind  speed  did  not  show  an  oscillating  pattern).  Figs.  C15  and  and  C16  show  that 
power  is  most  closely  correlated  vwth  Ow  (0.46)  and  power  is  negatively  correlated  with  both 
wind  speed  (-0.42)  and  vertical  wind  shear  (-0.26). 

In  summary,  these  cases  illustrate  that  while  there  is  a  relationship  between  wind  speed 
and  power  (Cn^)  on  a  large-scale  basis,  the  correlation  of  power  with  a*  (an  indicator  of  gravity 
wave  intensity)  is  much  more  robust.  The  vertical  wind  shears  associated  with  these  gravity 
waves  are  superposed  on  the  vertical  shears  of  the  prevailing  winds  and  lead  to  local  dynamic 
instabilities  that  break  down  into  turbulence.  This  scenario  follows  that  encapsulated  in  the 
model  of  Wamock  and  VanZandt  (1985)  and  is  similar  in  concept  to  that  given  in  Dewan  et  al. 
(1993),  and  follows  from  the  well-known  relationship  between  Cn^  and  Lq  plus  that  between  Lq 
and  micro-scale  wind  shear.  In  view  of  this  scenario,  it  is  clear  that  any  efforts  to  predict  C„ 
levels  must  take  into  account  the  strengths  of  gravity  wave  source  mechanisms,  i.e.,  it  is  not 
sufficient  to  predict  the  large-scale  flow  in  the  region  where  a  forecast  of  is  desired.  Wamock 
et  al.  (1988)  found  that  the  Cn^  model  developed  using  observations  fi’om  Colorado  (a  region  rich 
in  terrain-launched  gravity  waves)  had  to  be  modified  significantly  to  match  observations  taken 
over  the  flat  terrain  of  Illinois.  Since  gravity  waves  upwelling  from  the  planetary  boundary  layer, 
launched  by  flow  over  rough  terrain  or  by  convection,  are  often  the  most  important  source  of 
vertically  propagating  gravity  waves,  prediction  schemes  must  include  forecasts  of  boundary 
layer  winds  and  stability. 
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